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There are hints, both from cosmology and from neutrino oscillation experiments, that one or two sterile neutrinos with eV masses are 
favored. We consider the implications of combining data from short baseline neutrino experiments with cosmological observations. 
In particular we consider cosmological models where the equation of state of dark energy and the curvature are allowed to vary. We 
find that the consistency of results from short baseline experiments and cosmology is dependent on the cosmological data sets that are 
used. When the two are consistent, we find that the inclusion of information from short baseline experiments only mildly affects the 
preferred values of the cosmological model parameters. A universe with a cosmological constant type dark energy and zero curvature 
is still fully consistent with the data we use. We also report that a thermalized sterile neutrino in the ~ 10 eV range is fully consistent 
with a wide range of cosmological data as long as information on the matter power spectrum in not taken into account. 
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1. Introduction 

The many links between the very small and the very large are 
some of the most fascinating aspects of cosmology. To under- 
stand the Universe and its history, we must know the funda- 
mental building blocks of matter and their properties. Both par- 
ticle physicists and cosmologists benefit from the connections 
between their subjects. For example, cosmological observations 
have for a number of years provided the most stringent upper 
bounds on the neutrin o mass scale, and may one da y be in a po- 
sition to determine it (Lesgourgues & Pas tor! 1201 2b . One could 
also look at the situation in the opposite way and say that a bet- 
ter knowledge and understanding of the neutrino sector would 
help us to improve constraints on cosmological parameters. 

The number of neutrino species whose masses are below 
the GeV scale and which couple to the Z*^ boson, i.e., inter- 
act weakly, was determined to be 2.984 + 008 from LEP data 
(iNakamura & Particle Data Groupll2010l) . If there are more neu- 
trino types than the three we already know about, they must be 
very heavy or only couple to gravity, or both. Neutrinos that do 
not participate in the weak interaction are k nown as steri le . They 
appear in the so-called seesaw mechanism (IZubed2004 for gen- 
erating low neutrino masses, and are typically very heavy there, 
much heavier than the electroweak scale, in order to explain how 
low the masses of the ordinary, active neutrinos are. However, as 
long as it only interacts gravitationally there are no a priori con- 
straints on the mass of a putative sterile neutrino. 

The question of the existence of sterile neutrinos with 
masses in the eV range has attracted a lot of attention 



l20TTallbl: iKopp et alJ [20TTL [loTl iGiunti & Lavede^ I20T1I) . 
There are s everal experiment s hinting at their existence 



in the last few years (IDonini et al.i i2012t iGiunti & Lavederl 
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(Aguila r et al.ll200UlAgu ilar-Are v"alo et al ll20Tot [Mention etaH 
i201 It lAbdurashitov et al.l 120061) . However, the situation is far 
from clear , since other experiment s disfavour light sterile neu- 
trinos, e.g. lArmbruster et al.l(l2002l) : lAdamson et al1(l2010l) . 

On the other hand there have been indications from cos- 
mological probes that a value for the effective number of rel- 
ativistic particles higher than the standard value of A^eff = 
3.046 ( ManganoetaLl l2005t |Riemer-S0rensen et al.' '2013") is 
preferred, the so-ca lled Dark Radiation (Giusarmaet al. 20 ij; 
iHamann et al.ll2010l) . Although there are other explanations, one 
or two eV-mass sterile neutrinos could be the cause of this result. 

A number of papers have considered the issue of adding one 
or two l ight sterile neutrinos to the cosmological matter bud - 
get, e.g.. Archidiacono et al.' (l2012h : lArchidiacono et al.l(l2013h : 
Hamann et al. (2011, 2010). Most of these studies, however, 
have added sterile neutrinos to the ACDM model, where the cos- 
mological constant is the dark energy. The nature of the dark en- 
ergy is unknown and the subject of large observational projects 
like the future Euclid mission (Amendol a et al.ll2012l) . In this 
paper we will therefore allow the equation of state of the dark 
energy to be a parameter to be constrained simultaneously with 
other cosmological parameters and neutrino properties. 

In a previous paper (iKristiansen & Elgar0vll2Oin) two of us 
looked into the consequences for cosmology if the reactor data 
really demand the existence of one or two eV-mass sterile neu- 
trinos. The aim of this paper is to improve upon this analysis by 
including the full likelihood function from short base line (SBL) 
neutrino experiments. 
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The structure of this paper is as follows. In section 2 we 
summarize the theoretical background and describe our method 
of analysis. The cosmological data are treated briefly since our 
method of analysis is standard. Section 3 contains our results, 
and we summarize and conclude in section 4. 



2. Cosmological models with sterile neutrinos 

Sterile neutrinos with keV masses ha ve been of great interest as 
dark matter candidates in cosmology ('Kusenko'2009') . Recently 
it has been suggested that one or two sterile neutrinos with 
masses of a few eV lie behind so me puzzling fea t ures in neu- 
trino oscillation experiments, see lAbazaiian et al.l (1201 2h for a 
review. 

An important question when we turn to the cosmological 
implications of these two scenarios is whether these light, ster- 
ile neutrinos were thermalized in the early universe. We will 
assume that they were, since several studies (Hama nn et al.l 
[2010; Melchiorri et al. 2009; Kainulainen 1990) suggest fliat 
this was the case for the masses and mixing parameters we 
consider, prov ided that the initial lepton asymmetry was small 
dHannestad et al. 2012). This means that standard relation be- 
tween the sum of the neutrino masses and their contribution to 
the cosmic mass density parameter applies. 

Adding hght sterile neutrinos may cause some problems 
with BBN since the increased relativistic energy density results 
in a larger neutron-to-proton ratio, and l eads to an increased 
He/H mass fraction. Recen t analyses jMangano & Serpicol 
1201 It iPettini & Co"ok3 1201 2h found that BBN constrains the 
number of relativistic degrees of freedom to be A^eff < 4 at 95 
% confidence. So this means that the 3 active-i-1 sterile model 
we consider in this paper is just within the bounds. 



3. Neutrino oscillation constraints 

We considered the data of SBL neutrino oscillation exper- 
iments in the framework of 3 + 1 neutrino mixing, assum- 
ing that the mass 1114 of the fourth neutrino is much larger 
than those of the three ordinary neutrinos which generate the 
squared-mass differences responsible for solar and atmospher ic 
neutrino oscillations (see iGonzalez-Garcia & Maltonil dBosh V 
In this case, SBL oscillations depend only on the squared- 



mass difference Am^j ^ m^. We fitted the v„ 



and 



Ve ap pearance data of the LSND (Aguilaret al. 2001), 
KAR MEN lArmbruster et alj l2002h. NOMAD (.Astier et al., 



I2003h. MiniBooN E ( Aguil ar-Arevalo et"ani2012h and ICARUS 
dAntonello et alj I2OIZ) experiments, taking into account the 
results of Ve and Ve dis appearance da ta (Giunti et al. 2012) 

i including the reactor ([Mention et al.l |201 ih and Gallium 
Giunti & Lavedeij|201 1|) anomalies) and the co nstraints on 
and disappearance ( Archidiacono et aT]|2013l) . We calculated 



the Bayesian posterior probability of Am^j through the marginal- 
ization over the mixing angles. This probability is used as a SBL 
prior in the MCMC analysis of cosmological data. 

The posterior from the SBL analysis is shown in Figure [T] 
The SBL oscillation experiments are sensitive to the squared 
mass difference between two neutrino mass eigenstates, Arn^^. 
Since we assume the three active neutrino species to be mass- 
less, we can use the simple relation nty = -^Am^j . 



4. Cosmological parameter estimation 

For the cosmological parameter analysis we use a modified ver- 
sion of the publicly available Markov chain Monte Carlo code 
Fortran 90 code CosmoMC (Lewis & Bridle 2002). We have 
studied two different basic cosmological models, one minimal, 
flat ACDMmodel, and one model that is extended to allow for 
a non-zero curvature and dark energy equation of state, w, dif- 
ferent from -1. The latter model will be denoted wCDM+k. We 
have assumed the three active neutrino species to be massless. 
Thus, in the ACDM case we have varied the parameters {a>h, Wc, 
6, T, tis. In 10'*'As}, and for the wCDM +k case we have added 

and w as free parameters. We refer to Lewis & Bri dle (2OO2I1 
for exact parameter definitions. 

Then we have added one extra, massive, sterile neutrino 
species, v,, to the models, and included the mass of this ster- 
ile neutrino, niy (measured in eV/c^), as a free parameter For the 
sterile neutrino mass we analyze two cases. One where we use 
a flat prior on niy, and let only the cosmological data constrain 
its allowed mass range, and one where we impose a prior on niy 
from the analysis of short baseline reactor neutrino oscillation 
experiments. 

To include the SBL results in the cosmological analysis we 
simply add the log likelihood values from SBL and cosmological 
data for each step in the MCMC analysis. 

The backbone of the cosmological d ata we use are the9 
year CMB data from the WMAP satelUte dBennett etalJl2012l: 
Hinshaw et al. 2012), and the likelihood code provided with the 
data release. In addition we use s upernova type la data (SNla) 
from SDSS (iKessler et al.l l2009l) and matter power spectrum 
(P(k)) data from the SDSS DR7 LRG sample dReid et alJl2()Tol) 
and baryonic acou stic oscillation (BAO) data from WiggleZ 
(lBlakeetalJl201ll) and SDSS (Perciv al et al. 2010 ). In addi- 
tion we impose a prior on the Hubble parameter (HST) from 
iRiess et aL l (l20Tll) and on the baryon content from nucleosyn- 
thesis (BBN) (Slbh^ = 0.022 + 0.002). 

We have performed the analysis with the following combi- 
nations of data sets: 

- WMAP9: Only WMAP9 data. 

- Extended: WMAP9 + SNla + BAO + HST + BBN 

- Extended+Pk: As Extended, but with galaxy power spectrum 
P(k) data from SDSS instead of the SDSS BAO data. 



5. Results 

In Figure|2]we show ID marginalized parameter probability dis- 
tributions obtained for the ACDM+1 v, model for both WMAP9 
and Extended data, with and without the SBL prior on niy. 

For the sterile neutrino mass we notice that the distribution 
peaks at zero mass when we only take cosmological data into 
account. For WMAP9 data only the probability distribution of 
ifiy has a long tail of low but non-zero probability, leaving niy ba- 
sically unconstrained. We will discuss this tail more thoroughly 
later. For Extended data the 95% region for niy reaches up to 0.85 
eV. When including the SBL prior the niy distribution peaks at 
around 1 eV, as this region gives a reasonable fit both to the SBL 
and cosmological data sets. We see that the SBL prior is con- 
sistent with the cosmological analysis using both WMAP9 and 
Extended data. However, when we include matter power spec- 
trum data and use the Extended -i-Pk data, the cosmological limits 
on niy are tightened to niy< 0.46 eV (95% C.L.), which we con- 
sider to be incompatible with the SBL prior We conclude that 
the strong sterile neutrino mass constraints from the inclusion of 
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Fig. 1. Constraints on the mass of a single sterile neutrino from reactor SBL experiments. We use this likelihood distribution as a 
prior in our cosmological models. As we assume massless active neutrinos, the sterile neutrino mass, iriy, shown here, is simply 

taken as JAm^^ , where Am^j is the squared neutrino mass difference measured in the SBL experiments. 




Fig. 2. Marginalized probability distributions for cosmological parameters in a flat ACDM+lv, cosmology. Blue, dashed lines: No 
SBL-prior on niy. Red, solid lines: Including a SBL prior on niy. Thin lines: WMAP9 data. Thick lines: Extended data. In the 
niy panel the green dash-dotted line show the limit on niy for Extended+Pk data. In the other panels the black, dotted lines show 
distributions for a flat ACDM model without any sterile neutrinos, for comparison. 



P(k) data and the inferred limits from the SBL analysis cannot 
be valid at the same time. Therefore we do not present results 
with both Extended+Pk data and the SBL prior. 

Regarding the other cosmological parameters we notice that 
the inclusion of a sterile neutrino shifts several parameter values 
significantly compared to a standard ACDM universe without 
massive neutrinos. For example, the Hubble parameter. Ho, and 
dark energy and dark matter contents are altered, which leads to 
a shift in the age of the universe to lower values when the ster- 
ile neutrino is included. We see only minor shifts in parameter 
ranges when the adding SBL prior on m^is added. The SBL prior 



does however have a small efifect on the values of, for example, 
H() and Qa, which are drawn toward smaller values, and, as a 
result, the age of the universe is shifted towards higher values 
with the SBL prior included. Using the Extended data set and 
the SBL prior for the ACDM+v^ universe, we end up with an 
age estimate of the universe given by 13.10 + 0.08 Gyr, which 
can be compared to 13.68 + 0.08 Gyr when using the same cos- 
mological data sets for a ACDM model without sterile neutrinos. 

The reason for these shifts can be seen in Figure |3] Here we 
show two-dimensional confidence intervals (95% and 68%) for 
a few parameter combinations when using the flat ACDM+lvj 



3 



J. R. Kristiansen et al.: Cosmology with sterile neutrinos 



a 



0.14 



0.135 



P 0.13 



0.125 




0.74 
0.72 
0.7 



0.5 1 
m 



1.5 



0.68 
0.66 








0.5 1 
m 



1.5 





Fig. 3. 2D probability contours for selected parameters in the case of a flat ACDM+lvj model using Extended data. The contours 
indicate 95% and 68% confidence intervals. Blue lines show results without any SBL prior on niy, while the red contours indicate 
the results when including the SBL prior in the analysis. 



model and Extended data. We see that there are degeneracies 
between myand Qa and Hq, and as a consequence also the age of 
the universe. 

Next, we turn to the extended cosmological model, 
wCDM+^, where we allow the equation of state of dark energy, 
w, and the curvature of the universe, Q^^ to vary. In Figure |4] 
we show the probability distributions for the different parame- 
ters in this case. As these models are very poorly constrained by 
WMAP9 data only, we only show results from the Extended and 
Extended+Pk data sets. 

For the distribution of niy analyzed with the Extended data 
set we see that the distribution falls off in an almost Gaussian 
manner up to my~ 2 eV, and that it retains a small, non-zero 
probability for higher masses. We will discuss this high-mass 
tail further later on. The SBL prior is fully consistent with the 
extended data set, and adding this prior, constrains the neutrino 
mass to be between 0.96 eV and 1.40 eV at 95% C.L. Again, 
including matter power spectrum data, using the Extended+Pk 
data set we get limits on niy that are incompatible with the SBL 
prior (iriy] 0.55 eV (95% C.L.)). Therefore we do not show re- 
sults where the SBL prior is added to the Extended+Pk analysis. 

For the other parameters we find shifts of the order of 1 cr by 
adding a massive sterile neutrino, compared to a model without 
sterile neutrinos. Adding a sterile neutrino the preferred model 
is drawn towards less dark energy density, and more dark matter 
The distributions of D.^ and w shift by less than Icr for all the 
combinations of models and data sets shown here. This may be 
surprising, as we have found in a previous work that a sterile 
neutrino in the 1 eV range gives a preference for w < -1. In 
that analysis, however, we did not make use of BAO data. The 
inclusion of BAO data puts an additional strong handle on w. 



which makes it strongly constrained, regardless of the inclusion 
of additional massive neutrinos. 

Adding the SBL prior on the neutrino mass does not affect 
the distributions significantly. We do however see small shifts in 
parameters like the scalar spectral index, n,, where the signifi- 
cance of ris < 1 increases from 2. Icr to 3. Oct with the inclusion 
of an SBL prior This is important for inflationary models of the 
early universe. 

In the wCDM+k+Vs model we infer an age of the universe of 
12.89 + 0.22 Gyr when using the extended data set and the SBL 
prior 

In Figure |5] we show 2D distributions for a few parameter 
combinations in the case of a wCDM+k model analyzed with 
the Extended data set, with and without the SBL prior We notice 
that there is basically no correlation between the sterile neutrino 
mass and w, but that there are correlations between iriy and Wj, 
Oa and Q.^ . In particular we see that the anti-correlation between 
ris and niy show why SBL prior results in a stronger preference 
for n 5 < 1. 

Again, we notice the tail of an allowed region towards my> 2 
eV for the analysis using the Extended data set. In Figure |5] the 
95% region is indeed extending beyond my= 3 eV. In Figure|6]we 
show the marginalized probability distribution of niy for different 
combinations of models and data sets, without using the SBL 
prior For the Extended+Pk data, the probability for niy > 3 eV 
remains zero. However, for both the WMAP9 and Extended data 
sets, a new high-probability region appears in the range ~ 4 eV - 
15 eV. Indeed, for the ACDM+v^model analyzed with WMAP9 
data only, the preferred value of niy peaks at 13 eV. When the 
Extended data is used, the significance of this peak is lowered, 
and with the inclusion of matter power spectrum data the high- 
mass peak vanishes. 
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Fig. 4. Marginalized probability distributions for cosmological parameters in a wCDM+^+Vj cosmological model. Line colors 
corresponding to Figure |2] As this extended model is poorly constrained by WMAP9 data alone, we only show results from the 
Extended analysis. For the Extended+Pk cosmological data (green dash-dotted lines) we only show results without a SBL prior 



Does this peak hint towards the existence of a sterile neu- 
trino with a mass in the 10 eV range? Not really. In both the 
WMAP9 and Extended data sets the only direct handle we have 
on matter clustering comes from the CMB power spectrum. 
Neutrino masses affect the C MB power spectrum in several ways 
dLes gourgues & P astorll2006 ). but one of the most important ef- 
fects stems from the shift of the time of matter-radiation equality 
(teq). Neutrinos with masses less than ~ 2 eV will still be rela- 
tivistic at this time. Therefore, increasing the fraction of dark 
matter in the form of neutrinos, fy, will shift 1^^ to later times, 
which in turn will affect the structure formation. However, neu- 
trinos with masses > 2 eV will akeady be not-relativistic at f^,^, 
and will therefore not shift teq relative to a pure cold dark matter 
model. In this regard, neutrinos in the ~ 10 eV range, corre- 
sponding to the high-mass peaks in Figure |6] will therefore act 
similar to a cold dark matter component when it comes to the 
effect of teq on the CMB. When my~ 10 eV, fy will approach 1 
for fully thermalized sterile neutrinos (the exact value of course 
depending on the total dark matter content), meaning that almost 
all of the dark matter is in the form of sterile neutrinos. Then the 
normal upper bounds on the dark matter content of the universe 
will apply, explaining why the niy probability distributions falls 
of to zero for my~ 15 eV. 

When we include the matter power spectrum, using the 
Extended-i-Pk data, the situation changes. Here the large thermal 
velocities of the neutrinos give the neutrinos a free-streaming 



scale that suppresses structure growth at small scales. The sig- 
nature of a universe with ~ 1 will therefore look dramatically 
different from a universe with fy ~ 0. 

In Figure |7] we show the CMB and matter power spectra for 
cosmological models with different masses of the sterile neutri- 
nos. The other cosmological parameters are set to concordance 
values. The model with niy- 10 eV model therefore corresponds 
to a universe with a fy close to 1 . For the CMB power spectrum, 
the model with my= 10 eV actually has a larger resemblance to 
the standard ACDMmodel than the models with my= 1 eV or 
;«v= 3 eV. This can be understood by the reasoning above, re- 
garding the neutrino mass effect on t^q. However, when we turn 
to the matter power spectrum, we see that the niy- 10 eV model 
gives an extreme suppression of small scale structure, in total 
disagreement with all matter power spectrum data. Also, such a 
high value of niy is in disagreement with the SBL results used 
here. 

Thus, we regard the high-mass peaks of niy not to be com- 
patible with our observed universe, at least not within our model 
framework. In the analysis of the data from WMAP9 and the 
Extended data sets shown above, we have therefore imposed a 
cutoff at niy = 5 eV in the analysis of the Monte Carlo chains, 
to avoid the inclusion of the region of parameter space with an 
artificially high niy. 



5 



J. R. Kristiansen et al.: Cosmology with sterile neutrinos 



1.02 
1 

0.98 
0.96 
0.94 








0.74 
0.72 
0.7 
0.68 
0.66 








m 



m 



a 




-0.9 
-1 
-1.1 
-1.2 
-1.3 








m 



Fig. 5. 2D 95% and 68% probability contours for selected parameters in a wCDM-i-A:+VjCOsmological model, using Extended data. 
Blue contours: No SBL prior. Red lines: With the SBL prior included. 




Fig. 6. Probability distributions for niy for different combinations of data sets and cosmological models without SBL prior. Red, 
solid line: ACDM+v., model, WMAP9 data. Black, dotted line: ACDM+v, model. Extended data. Blue, dashed line: wCDM-i-A; 
model. Extended data. Green, dash-dotted line: wCDM+^ model, Extended-t-Pk data. 



6. Discussion and conclusions 

Our analyses show consistency between cosmology and SBL 
data when the cosmological data consists of the combination 
WMAP9-SNla+BA0+HST+BBN. The small peak in the SBL 
prior seen around 2.5 eV is removed by cosmology, but the 
higher peak around 1 eV remains. However, cosmology becomes 
inconsistent with SBL when the BAO data are replaced with 
the galaxy power spectrum estimated from the SDSS LRG sam- 
ple. The latter result depends on the relation between the dis- 
tributions of galaxies and the dark matter, a relation which is 
non-trivial and not understood in detail. Whether a better under- 
standing of the small-scale clustering of galaxies and matter will 



leave room for a sterile neutrino with mass ~ 1 eV remains to be 
seen. There is, however, no reason to doubt that the galaxy power 
spectrum data exclude the peak around 10 eV in the probability 
distribution for niy that appears when using CMB + background 
evolution data only. 

While iKristiansen & Elgar0vl ( 1201 Ih found a preference for 
the equation of state of dark energy w < - 1 when a sterile neu- 
trino was included, in this work we find that the cosmological 
constant lies within Icr contour. The BAO data play a crucial 
role for this result by effectively removing the degeneracy be- 
tween My and w. 

Our results are consistent with those of (lArchidiacono et all 
i2013i) . They considered a wider range of cases for sterile neu- 
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Fig. 7. CMB and matter power spectrum for ACDM and ACDM+v, models with different values of niy. The total dark matter density 
is kept constant, so the plots show the effect of changing different amounts of cold dark matter into hot dark matter in form of fully 
thermalized sterile neutrinos. The other cosmological parameters are kept to standard, concordance values. Red, solid line: ACDM 
model. Black, dotted line: ACDM+v, model with niy-l eV. Blue, dashed line: ACDM+v, model with my=3 eV. Green, dash-dotted 
line: ACDM+Vjmodel with niy-lO eV. 



trinos, e.g., both 3+1 and 3+2, and the possibility of incomplete 
thermalization, but did not allow the equation of state of dark 
energy to vary, and assumed a spatially flat universe. We have 
allowed w and the spatial curvature to vary, but have considered 
the 3+1 case with a fully thermalized sterile neutrino. Both anal- 
yses agree that a sterile neutrino with mass ~ 1 eV is consistent 
with cosmology. 

The CMB on its own ha s a limited ability to constrain neu- 
trino masses dLesgourgues & Pastor 2006), even with the pre- 
cision we can expect from the soon-to-be released Planck data. 
This is well illustrated by the peak around 12 eV for the sterile 
neutrino mass we obtained when not using data on the matter 
power spectrum. But by removing parameter degeneracies when 
combined with other probes, we can expect better constraints on 
niy. It will be interesting to see if cosmology can be consistent 
with a 1 eV sterile neutrino when Planck data are combined with 
a cleaner probe of the clustering of matter than the galaxy power 
spectrum. 
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